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We study the magnetization plateau at 1/4 of the saturation magnetization of the S = 1 
antiferromagnetic spin ladder both analytically and numerically, with the aim of explaining 
recent experimental results on BIP-TENO by Goto et al. We propose two mechanisms for the 
plateau formation and clarify the plateau phase diagram on the plane of the coupling constants 
between spins. 
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The magnetization plateaux of low dimensional spin 
systems have been attracting much attention be- 
cause of its full quantum nature. Very recently, 
a new organic tetraradical, 3,3',5,5'-tetrakis(Af-ieri- 
butylaminoxyl)biuhenyl, abbreviated as BIP-TENO, has 
been synthesized. Ef This material can be regarded as an 
5 = 1 antiferromagnetic two-leg spin ladder, as explained 
later. Goto et al& measured the magnetization curve 
of BIP-TENO at low temperatures in pulsed high mag- 
netic fields up to about 50 T. They found M ~ up to 
about H c i ~ 10 T, suggesting the non-magnetic ground 
state, which was consistent with the result of suscepti- 
bility measurement. Another remarkable nature of their 
magnetization curve is the plateau at M = M s /A above 
H C 2 ~ 45 T, where M s is the saturation magnetization. 
Unfortunately the end of the M s /4 plateau is unknown 
because of the limitation of the strength of their mag- 
netic field. 

Sonsidering the structure of BIP-TENO, Katoh et 
proposed its model of Fig. 1. Since the ferromagnetic 
coupling Jp is thought to be much larger than other anti- 
ferromagnetic couplings J' AF and Jafi, two spins coupled 
by Jp effectively form an S = 1 spin. Thus the model of 
Fig. 1(b) is reduced to an S = 1 antiferromagnetic two- 
leg ladder, as shown in Fig. 1(c). The relation between 
the coupling constants of these two models are 
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In the following, we consider the model of Fig. 1(c) to in- 
vestigate the M s /4 magnetization plateau of this model 
at zero temperature. The plateau at M = is inter- 
preted as the manifestation of the Haldane state. On the 
other hand, the M s /4 plateau is not easily understood. 
The necessary condition for the quantized magnetization 
plateau by Oshikawa, Yamanaka and AffleckcP is 



Q(S U 



(m)) = integer, 



(2) 



where Q is the spatial period of the wave function of the 
state measured by the unit cell, S'unit and (m) are the 
total spin and the magnetization per unit cell, respec- 




Fig. 1. (a) Schematic illustration of BIP-TENO. Open circles 
represent S = 1/2 spins, (b) Plane type illustration of BIP- 
TENO. (c) Effective model of BIP-TENO in case of J F -» oo. 
Open ellipses represent 5 = 1 spins. 



tively. Applying the theorem to our model, S'unit = 2 and 
m = 1/2, this condition is not satisfied if Q = 1. Thus, 
Q = 2 is required at least, which means the occurrence of 
the spontaneous symmetry breaking in the M s /4 plateau 
state. The magnetization plateaux caused by the spon- 
taneous symmetry breaking were discussed for the M s /2 
plateau of the 5=1/2 zig-zag ladderoBEl' (equivalent to 
the bond-alternating chain with next-nearest-neighbor 
intej^ticms), the M s /2 plateau of the 5 = 1/2 lad- 
derJiJ&BEl" and also the ftA)M s plateau of the 5=1/2 
distorted diamond chainJid'LL 2 ? and so on. 
The Hamiltonian of the model of Fig. 1(c) is 



1=1,2 j=l 



3 = 1 



2 



Author Name 



(3) 



1=1,2 j=l 



where S denotes the spin-1 operator, j the rung number 
and 1 = 1,2 the leg number. The last term is the Zeeman 
energy in the magnetic held H. When Jj_ = 0, our model 
is reduced to that of two independent usual S = 1 chains, 
in which no plateau is expected in the magnetization 
curve except at M = 0. 

Let us consider the opposite limit J± 3> J\ by use of 
the degenerate perturbation theory. When Jj = 0, all 
the rung spin pairs are mutually independent. In this 
case, at M s /A, half of the rung spin pairs are in the state 

V(0,0) = (l/v / 3)(|U> + ||T>-|00)) (4) 

and the remaining half pairs are in the state 

^(1,1) = (1/V2)(|T0) + |0T)) (5) 

where ip(Stot, S^ ot ) is the wave function of two S = 1 
spins coupled by antiferromagnetic interaction J± with 
the quantum numbers S'tot and S* ot . These wave func- 
tions have lowest energies in the subspace of S* ot = 
and S* ot = 1, respectively. Other possible selections have 
higher energies. The M s /4 state is highly degenerate as 
far as Ji = 0, because there is no restriction for the con- 
figurations of these two states. This degeneracy is lifted 
up by the introduction of J\ . To investigate the effect of 
Ji, we introduce the pseudo-spin T with T = 1/2. The 
| ff) and | -IJ.) states of the T spin correspond to ^(1,1) 
and ip(0,0), respectively. Neglecting other seven states 
for the rung spin pairs, the effective Hamiltonian can be 
written as 
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in the lowest order of Ji , where 



J cff 



(8/3)Ji 



J. 
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(6) 



(7) 



Thus the M s /4 plateau problem of the original model is 
mapped onto the M = problem of the usual T = 1/2 
antiferromagnetic XXZ spin chain with the nearest- 
neighbor (NN) interaction. As is well known, depending 
on whether A G ff = J<f ff /J^ < 1 or A G ff > 1, the ground 
state of eq.(||) is either the spin- fluid state (gapless) or 
the Neel state (gapful), which correspond to the no- 
plateau state or the plateau state of the original model, 
respectively. 

We see A cff = 3/16 from eq.(0), which means no M s /4 
plateau in the model of Fig. 1(c), at least when J± ^> 
J\. As already stated, there is no M s /4 plateau also 
when Jl = 0. Thus the M s /4 plateau will not appear 
in the whole region of Ji/Ji for the model of Fig. 1(c), 
although we cannot completely exclude the possibility of 
the M s /4 plateau at the intermediate region. We note 
that the absence of the M s /4 plateau is also established 
by numerical method, as shown later in Fig. 4. 



:s for the M s /2 plateau of the 5=1/2 
i " it has shown that the second and/or 
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Fig. 2. 5 = 1 ladder model with second and third neighbor in- 
teractions. 



Considering this fact, let us introduce the second (J2) 
and third (J3) neighbor interactions to our model, as 
shown in Fig. 2. Including these interactions, the effec- 
tive Hamiltonian (|^) is modified to the generalized ver- 
sion of the T =1/2 XXZ Hamiltonian with the NN and 
the next-nearest-neighbor (NNN) interactions: 

^eff - £ {J%(T?1f +1 + 7JJJ +1 ) + J* S T*T* +1 } 
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where 



rj = (8/3)07, 



J- 



ff — v°/ "A"l ~ ^2), 
f 2) = (8/3) J 3 , 



= (1/2)07! + J 2 ) (9) 



J e f } = (1/2) J 4 



(10) 



Note that the XXZ anisotropy parameters of the NN 
and NNN interactions are different with each other jn 
general. In this model, three phases are expected :li3'tLS ) 
these are the spin-fluid state (gapless), the Neel state 
(gapful) and the dimer state (gapful). 

For simplicity, first we consider the J2 > 0, J3 = case. 
In this case, the NNN interaction in eq.(8) vanishes and 
the anisotropy parameter is 

3 Ji + J 2 ri1 > 

A eff = 7- (11) 

16 Ji - J 2 

Therefore the Neel condition A er f > 1 is satisfied when 
J2/J1 > 13/19. This type of M s /4 plateau is due to the 
Neel mechanism and called plateau A from now on. 

Next we consider the J2 = 0, J3 > case. In this case, 
the anisotropy parameters of the NN and NNN interac- 
tions have the same value A Gff = 3/16 and J^ (2) /J^ = 



7-2(2) , „ 

•^eff / J cff 



J3/J1 



From the phase diagram of S 
XXZ chain with the NN and NNN interactions! 
the ground state of (8) at A c g = 3/16 is either the 
spin-fluid state or the dimer state depending on whether 
J3/J1 < 0.31 or J 3 /Ji > 0.31. This type of M s /4 
plateau, which is brought about by the dimer mecha- 
nism, is called plateau B. 

The physical pictures of these two plateau states are 
shown in Fig. 3. As is clear from Fig. 3, both plateau 
states are doubly degenerate due to the spontaneous 
symmetry breaking. Then the period of the state is 
Q = 2 in both plateau states, from which we see the 
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Fig. 3. (a) Physical picture of the plateau A (Neel mechanism). 
Open and shaded ellipses represent i/>(0, 0) and V(b 1) states of 
the rung pairs (| ij-) and | ft) in the T picture), respectively, (b) 
Physical picture of the plateau B (dimer mechanism). 



realization of the necessary condition eq.(g). The tran- 
sition from the no-plateau state to the plateau A or B 
state (in the language of T, from the spin-fluid state to 
the Neel jOTpdimer state) is.af.the Berezinskii-Kosterlitz- 
ThoulesskJy> (BKT) type.B0J In general case J 2 , J 3 > 
0, the transition from the plateau A state to the plateau 
B state, may occur, which has the Gaussian universal- 
ity.BB> 

We have performed numerical diagonalization by 
Lanczos method up to 16 spins for the S = 1 model with 
Ji and/or J 3 to confirm to existence of two plateaux, pre- 
dicted by use of the degenerate perturbation theory. As 
already stated, the transition from the no-plateau state 
to the plateau A or B state is of the BKT type, of which 
boundary is very difficult to determine from the numer- 
ical data when the conventional methods are used. This 
is mainly due to the logarithmic size corrections associ- 
ated with the BKT transition, We have used the level 
spectroscopy methodjSBIiJ'EZloES' by use of which the 
BKT phase boundary can be determined free from the 
most dominant logarithmic size corrections. In the level 
spectroscopy method, we use three excitations defined by 
(for simplicity we are writing the L — An case, n being 
an integer) 
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(12) 
(13) 
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where Eo(L,M,k) means the lowest eigenvalue of the 
S = 1 ladder Hamiltonian in the subspace where the 
eigenvalue of the operator 5 t 2 ot ee J^V Sfj is M with 
the momentum k and system size L = An under the pe- 
riodic boundary condition. En{L, L/2,0) is nothing but 
the ground state energy. Eqx is the lowest energy with 



M = L/2, k = 7r subspace having the eigenvalue P = — 1 
for the space inversion operation of the rung number 
i —> L — i + 1. This excitation is called Neel excitation. 
We have P = — 1 for Eob, which is named dimer excita- 
tion. In the no-plateau state, plateau A state and plateau 
B state, the lowest excitations should be Ai, AOA and 
A b, respectively. Thus the boundaries between these 
three phases can be determined from the crossing points 
among these three excitations with sweeping the coupling 
constants.E-f 
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Fig. 4. The Ms/4 plateau phase diagram (a) on the Ji — J2 plane 
when J3 = 0, and (b) on the J\ — J3 plane when J2 = 0. Here 
Ji = Ji/Jj_ for i = 1,2,3. The closed diamonds represent the 
BKT phase boundaries and open circles the first order phase 
boundaries. 



By use of the level spectroscopy, we obtain the M s /A 
plateau phase diagrams in Fig. 4. The estimated error 
bars of the boundary points are much smaller than the 
size of marks. The line symmetry with respect to the 
line Ji = Ji in Fig. 4(a) is due to the equivalence of the 
model for interchanging J\ <-> J2 when J3 = 0. The slope 
of the BKT lines near the origin are very close to the 
predicted values, 13/19 in Fig.4(a) and 0.31 in Fig.4(b), 
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Fig. 5. (a) Classical picture of (a) the no-plateau 1 state, and (b) 
no plateau 2 state. 



which confirms the high reliability of our analysis of the 
numerical data by use of the level spcctroscopyjnciJiod. 
In Fig. 4(a), there are two no-plateau phases ptrH of 
which classical pictures are shown in Fig. 5. Note that 
the long range order is destroyed by quantum fluctua- 
tions in no-plateau states. The boundaries between two 
no-plateau states are determined by the crossing of the 
ground state. 

Following the present study, the critical value of J3 
for the realization of the plateau when J2 = is smaller 
than that of Ji when J3 = 0, irrespective of the value of 
Ji . Furthermore the plateau B phase can appear for any 
value of Ji, although the plateau A cannot be realized 
for larger values of J\. Thus the plateau B is easier to 
be realized than the plateau A. 

One may think that there should exist dashed-line in- 
teractions in Fig. 6(a). The model of Fig. 6(a) has four 
bond-alternating S = 1/2 chains. At a glance, the unit 
cell seems to consist of eight 5 = 1/2 spins, which may 
brings about the M s /4 plateau without any spontaneous 
symmetry breaking. However, the model of Fig. 6(a) can 
be re-drawn in the plane form of Fig. 6(b), indicating 
that the unit cell is still composed of four S = 1/2 spins. 
Thus this new interaction cannot be the direct reason for 
the M s /4 plateau. In fact, in the degenerate perturba- 
tion theory used in this letter, this new interaction only 
modifies Jfg and J^ ff in eq.(7) without changing their 
ratio A e ff = J^/J* s . A similar situation can be found 
in two-leg boncLalternating ladder investigated by Cabra 
and Grynberg .Ef 

In conclusion, we have investigated the M s /4 plateau 
of the S = 1 spin ladder, proposing two mechanisms 
for the plateau. We have also clarified the plateau-no 
plateau phase diagram by use of the degenerate pertur- 
bation theory, and the level spectroscopy analysis of the 
numerical diagonalization data. The full study of the 
magnetization process of the present model will be pub- 
lished elsewhere. 
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Fig. 6. (a) Picture of four bond-alternating 5 = 1/2 chains, (b) 
Plane type illustration of (a). The unit cell consists of four S = 
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